Interacting Anyonic Fermions in a Two-Body Color Code Model 
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We introduce a two-body quantum Hamiltonian model of spin- 1/2 on a 2D spatial lattice with 
exact topological degeneracy in all coupling regimes. There exists a gapped phase in which the 
low-energy sector reproduces an effective color code model. High energy excitations fall into three 
families of anyonic fermions that turn out to be strongly interacting. The model exhibits a Z2 x Z2 
gauge group symmetry and string-net integrals of motion, which are related to the existence of 
topological charges that are invisible to moving high-energy fermions. 
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I. INTRODUCTION 

The Kitaev model on the honeycomb lattice [l[ has 
attracted a great deal of attention in the recent years 
0, H, 0, IE] since it offers the opportunity to study many 
properties of topologically ordered systems in a very well- 
suited scenario for condensed matter. One of the most 
interesting properties of the Kitaev model is that one of 
its phases effectively reproduces the famous toric code. 
This is the first example of topological quantum error 
correction and plays a major role in quantum informa- 
tion. 

Yet, there is another family of topological codes, color 
codes, which exhibit very remarkable properties. They 
allow an appropriate implementation of the whole Clif- 
ford group of quantum gates, which are essential in quan- 
tum information tasksjfj. Topological color codes (TCC) 
are constructed with quantum lattice Hamiltonians that 
typically demand 6-body terms in a 2D spatial lattice, re- 
quirements quite unrealistic in a condensed matter frame- 
work. Up to now, it has remained a challenge to find a 
2-body Hamiltonian capable of hosting a TCC in a par- 
ticular coupling regime. 

In this paper we provide such a 2-body spin-i Hamil- 
tonian model. It turns out to exhibit very rich physics 
and, quite remarkably, it does not belong to the family 
of models originated after Kitaev's model. These lat- 
ter models are defined on trivalent lattices, which allows 
a fermionization yielding a free fermion exact solution. 
Instead, the lattice of our model is 4-valent, a sharp dif- 
ference that prevents complete solvability and gives rise 
to very interesting features not present in the mentioned 
models: 

i/ Exact topological degeneracy in all coupling regimes, 

rooted on the existence of string-net integrals of motion. 

This degeneracy, related to certain 'invisible' topological 

charges, is 4 9 -fold in surfaces of genus g. 

ii/ Emergence of three families of strongly interacting 

fermions with semionic relative statistics. 

iii/ An exact Z2 x Z2 gauge symmetry. Each family of 

emergent fermions sees a different Z2 gauge subgroup. 

Thus, the model admits an exact analysis of many inter- 



esting properties. This is so mainly because there exist 
local integrals of motion in an amount of w of the total 
number of spins. 

Although we will focus on a particular phase of the 
model, namely the one that effectively yields TCCs, a 
rich phase diagram beyond this gapped phase is to be 
expected in analogy with [l[ . This includes the possibility 
of non-Abelian anyons and other phases with interesting 
many-body effects. In addition, it is possible to break a 
symmetry of the model, which we call color symmetry, 
whereas the exact features above are kept. This paves 
the way to a yet more complex phase diagram. 

The properties of the model make it a good candidate 
for an experimental realization, by means of some engi- 
neering scheme like polar molecules on optical lattices [8] . 
It is also amenable to numerical computations Q with 
several methods, which would help to uncover some of 
its non-perturbative aspects and phases. 

In summary, the model represents a relevant contri- 
bution to the difficult task of searching for systems with 
emerging anyons and topological order, which are concep- 
tually important phenomena and a source of new physics. 

This paper is organized as follows: in Sect. |TT]we intro- 
duce the quantum Hamiltonian model based solely on 2- 
body interactions between spin-i particles. The lattice is 
two-dimensional and has coordination number 4, instead 
of the usual 3 for the Kitaev model. It is pictured in Fig.Q] 
and it is called ruby lattice. In Sect. IIII[ we introduce 
a mapping from the original spin-i degrees of freedom 
onto bosonic degrees of freedom in the form of hard-core 
bosons which also carry a pseudospin. In Sect. IIVI we 
describe the constants of motion of our model, specifying 
both its various types (normal strings and non-standard 
stringnets), as well as its number of | of the total num- 
ber of spins. With the help of the previous mapping, 
it is possible to analyse very interesting qualitative fea- 
tures of our model at non-perturbative level. In Sect. [Vj 
we show how the emerging quasiparticle excitations of 
our model are anyonic fermions with strong interactions, 
which are a manifestation of the string-net structure of 
the constants of motion. Thus, the model is not a free 
fermion model. In Sect. IVII we solve the problem of 
finding a 2-body quantum lattice Hamiltonian hosting 
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FIG. 1: Three different points of view of the system, (a) is the 
physical one, with vertices representing spins and colored links 
representing a^a^, a^aj^ and o z a &l two-body interactions. In 

(b) both spins and hexagonal plaquettes have been colored. 

(c) is the reduced lattice A, where vertices represent sites, 
that is, the blue triangles of (a) . We show a magnified site in 
(a,b) and a closed string in (c). 

the topological color code as one of its phases. Interest- 
ingly enough, the relevant properties of the TCC remain 
valid at the non-perturbative level as well. We also de- 
scribe the notion of invisible charges and discuss anyon 
condensation, and its implications. Sect. IVIIl is devoted 
to conclusions. 



II. HAMILTONIAN MODEL WITH 2-BODY 
INTERACTIONS 

The model of interest lives in the lattice of Fig. [Ha). 
Notice that, unlike in p], 0, 0, H|, it has coordination 
four. Vertices represent spin-i systems and links two- 
body interactions. The Hamiltonian is 

{x, red links 
y, green links (1) 
z, blue links 

where the couplings J w are real, a w are the Pauli 
matrices and the sum extends over all links. Notice that 
blue links are special since they form triangles, which we 
will refer to as sites for reasons that will be apparent 
below. 

Sites are the vertices of a reduced lattice A in which 
links are given by the pairs of parallel red and green links 
of the original lattice, see Fig. [Tfc). This reduced lattice 
will play an essential role in the study of the properties 
of the model. It is an hexagonal lattice and thus has 3- 
colorable plaquettes. We choose to color them with red 
(r), green (g) and blue (b). We also color links accord- 
ingly, in such a way that c-links connect c-plaquettes. 
Here and in what follows, we use the letter c to denote 
color variables. We remark that this coloring has noth- 
ing to do with the one used for links in Fig. QJa). The 
coloring of the reduced lattice induces a coloring of some 
plaquettes of the original lattice, see Fig. QJb), and also 
a vertex coloring that labels the spins within a site. 



Although we have given a particular lattice for con- 
creteness, the analysis that follows is much more general. 
Instead of having a hexagonal reduced lattice, it could 
be any trivalent lattice with 3-colorable plaquettes. Such 
lattices were named 2-colexes in Any closed sur- 

face will work, orientable or not. That the construction 
works on non-orientable surface is a consequence of the 
equivalence of vertex bicolorability and orientability in 
2-colexes Q. 

III. BOSONIC MAPPING 

In section IVII we will analyze in detail the regime 
J z > 0, | J x \, \ J y \ <C J z . Let us set for simplicity J z = 1/4. 
Then in the extreme case J x — J y = the system consists 
of isolated triangles, one per site. In an energy eigenstate, 
each of them contributes an energy —3/4 or 1/4, so that 
we can attach a quasiparticle with energy gap equal to 
1 to each triangle. With this motivation, we give here a 
map to a new system in which these quasiparticles are 
explicit. The mapping is exact, so that only the physi- 
cal picture is changed. It has the advantage of isolating 
those degrees of freedom that survive once the integrals 
of motion to be described below have been fixed. 

In Fig. [1] (b) each spin in a site has been identified 
with a color. We label the corresponding Pauli operators 
as a™ with c = r,g,b. Consider the Hilbert space He 
with orthonormal basis {|0), |r), |g), |b)} and introduce 
the colored annihilation operators 

6 c :=|0)(c|, c = r,g,b, (2) 

so that TLq represents a hardcore boson with three possi- 
ble color states. The number operator n and the colored 
number operator n c are 

n := ^ n c , n c :=b\b c . (3) 

c 

At each site, we attach such a boson and also an effective 
spin-i. We have to relate this degrees of freedom to the 
original three spins in the site, which are colored as in 
Fig. [ljb). The mapping can be expressed by relating 
bases of both systems. In particular, taking the usual 
up/down basis for the three physical spins and the tensor 
product basis 

\a,d) = |a}® \d), a=t,|, rf = 0,r,g,b, (4) 

for the effective spin and hardcore boson system we have 



T,o) = 


1 TTT), 


l,o}= 


1 Hi), 


(5) 


T,r) = 


1 Hi), 


U,r> = 


1 ITT), 


(6) 


UN 


I in), 


1 l,s) = 


1 TIT), 


(7) 


t,b) = 


I lit), 


U,b)= 


I m>. 


(8) 



More compactly, the mapping can be expressed by iden- 
tifying operators as follows 

a z c = t z ® p c , a v c =t v <g> (b\ + b c + s v r c ), (9) 
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where v = x, y, s x :— —s y := 1, the symbols r denote 
the Pauli operators on the effective spin and we are us- 
ing the color parity operators p c and the color switching 
operators r c defined as 

Vc ■= 1 - 2(n s + ng), r c := btbs + blb s , (10) 

where the bar operator transforms colors cyclically as 
follows: 

f:=g, g := b, b:=r. (11) 

From this point on we will be working always in the 
reduced lattice A. For compactness, we will use a sim- 
plified notation in which site indices are supressed and 
only relative positions are indicated: the notation O jC 
means O applied at the site that is connected to a site of 
reference by a c-link. For example, if i,j,k are neighbor- 
ing sites, with i the reference site and j, k are connected 
to i by a g-link and a r-link respectively, then instead 
of AiBjCk we write AB g C r . That is, we only indicate 
the relationship between sites, and completely omit the 
reference site i, which is therefore implicit. Also, we in- 
dicate the w — x,y,z indices of J w , t w , s w in terms of 
two colors as follows: 

c|c := z, c\c :— x, c\c :— y. (12) 

Then, the Hamiltonian {1} can be exactly transformed 
into a new form 

H = -37V/4 + Q - ]T ]T J c > lc T L c \ (13) 

A c^c' 

with N the number of sites, Q :— n the total number 
of hardcore bosons, the first sum running over the N sites 
of the reduced lattice, the second sum running over the 
6 combinations of different colors c, d and 

/ f c ' 4- v c> r c ' 

Tl =< + ^-^ + ^ + h.c, (14) 

a sum of several terms for an implicit reference site, ac- 
cording to the notation convention we are using. The 
meaning of the different terms in (|14|) is the following. 
The operator t c c is a c-boson hopping, r c c switches the 
color of two c- or c-bosons, u c c fuses a c-boson with a c- 
boson (or a c-boson) to give a c-boson (c-boson) and v c c 
destroys a pair of c-bosons. A pictorial representation is 
offered in Fig. Ufa). The explicit expressions are 

<-=4b c b\ c „ ri :=T*'r c r CtC ,, (15) 

U c c := s C '| c t c c 6 c r CjC /, v c c := r c c b c b c . d , (16) 

where we are using the notation 





FIG. 2: Several c-fermion processes, in the effective hexagonal 
lattice, c-fermions are depicted as c-colored balls, and their 
movement with c-colored directed strings. When the initial 
and final color charge at a site are different before and after 
a process, we show the transition with an small black arrow, 
(a) Examples of the effects of the terms in (|13|l . From top to 
bottom, the corresponding terms are t T s , itj 3 , r t an d Wg. (b) A 
typical c-fermion interaction process, with hoppings, fusions 
and splittings, (c) The exchange of two identical b-fermions. 
(d) An r-fermion surrounds a region A. The phase that it 
picks up depends both on the number of g- and b-fermions 
in A and the vortex states of the plaquettes marked with an 
spiral. 

IV. CONSTANTS OF MOTION 

In this section we explore several integrals of motion 
of the Hamiltonian (fT3")) . 

A. Color charge 

A basic property of the Hamiltonian fT5|) is the exis- 
tence of a Z2 x Z2 charge. Let us label the elements of 
Z2 x Z2 with colors as {e, r, g, b}, where e is the identity 
element. We also label its irreps as {x e , Xr, X g , Xb} set- 
ting Xe(c) := Xc(c) := ~Xc(c) := 1. Let us arbitrarily 
attach a charge to each family of hardcore bosons. In 
particular, we attach the irrep Xc to each c-boson. Let 
Q c = J2a nc b e t ne total number of c-bosons. The total 
color charge 

XA := X? r X^X? b (18) 

is preserved by the Hamiltonian. This can be checked 
directly or noting that the operators 

Xa(c) = =JJPc (19) 

A 

commute with all the terms in IT4]) . In section [V] the 
meaning of color charge will become clear, when we relate 
it to a Z2 x Z2 gauge field. 

B. Plaquette operators 

We describe now certain plaquette operators that com- 
mute with each other and with each of the terms of the 
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Hamiltonian. As we will show below, the corresponding 
degrees of freedom should be regarded as vortices. For 
each color c and each plaquette tt there is a constant of 
motion of the form 

s%:= n^ ic Vc<* c , (20) 

where d is the color of the plaquette tt, the product runs 
through its sites and * is just a convenient symmetric 
color operator defined by 

c*c:=c, c-k c := c* c :— c. (21) 

This operators are not all independent. They are subject 
only to the following constraints 

n^=(-!) s/2 . n ^=(-!) s/2 . ( 22 ) 

7rGA c— r,g,b 

where s is the number of sites of a given plaquette tt, and 
we are supposing in the first equation that the lattice 
forms a closed surface. It follows that there are 2F — 2 
independent plaquette operators, with F the number of 
faces, or plaquettes, of the reduced lattice A. The total 
color charge is not independent of plaquette operators, 
because 

Xa(c) = J] S E n I] Sl (23) 

where A c denotes the subset of c-plaquettes of A. 

In section |V] we show that plaquette degrees of free- 
dom can be regarded as vortices, as they correspond to 
a Z2 x Z2 gauge field, which has therefore no dynamics. 
The correspondence between plaquette operator eigen- 
states and group values is as follows. First, for each 
c-plaquette we introduce the following alternative pla- 



quette operators: 




K 


JX ^7T' 


K 


._ oS/2 OE 
Jy j tt: 


K 




where j w := J W /\J W \ 


. The element j, 6 Z2 x Z2 



we attach to each plaquete tt is determined, for a given 
eigenstate of the operators (|2"4"|) . by the conditions 

XM = B%. (25) 

These equation always have a solution because (B^) 2 = 
B^B&B^ = 1. Given a region or collection of plaquettes 
A, we will use the notation g\ := JlvreA 9™- 

C. String operators 

Plaquette constants of motion can be generalized to 
closed strings and, for that matter, also to string- nets 0- 



For a closed string we mean a connected path in the 
lattice with no endpoints, see Fig. [He). To any such 
string 7 and color c we attach a string operator 

^:=nr c|c W, (26) 

7 

where the product runs through the sites of 7 and, for 
each site, d is the color of the plaquette arround which 7 
turns at the site. This is exemplified in Fig. []Jc), where 
we have marked with a dot the corresponding plaquette 
for each vertex of the string. As in the case of plaquette 
operators, we have (S^) 2 = (— 1)' S!yS c S^ — 1 with s the 
number of sites of 7. Closed string operators commute 
with all plaquette operators and Hamiltonian terms, but 
not always with each other. Namely, if the strings 7, 7' 
cross once then [S°,S°,] = but {S°, S°,} = 0. 

This anticommutation property, which remarkably is 
not present in turns out to be crucial. In systems 
with nontrivial topology it is the source of an exact de- 
generacy of the Hamiltonian. For example, the lattices of 
Fig. [1] with periodic boundary conditions live on a torus, 
which has two independent nontrivial loops 7, 7'. Then 
anticommutes with Sy and commutes with and 
SL, showing[ll| that the Hamiltonian is at least 4-fold 
degenerate. More generally, in a closed surface with Euler 
characteristic x the degeneracy is 2 2_x . In particular, in 
orientable surfaces of genus g there exist a 4 9 -fold topo- 
logical degeneracy. This result is exact and independent 
of the particular phase of the system we are in. In or- 
der to label the global flux degrees of freedom attached 
to nontrivial string operators we can use the eigenvalues 
of 2 — x topologically nontrivial and independent closed 
string operators of a given color. 

V. ANYONIC FERMIONS 

In the previous section we have learned that the sys- 
tem is divided into sectors with a given vortex and global 
flux configuration. We now want to understand the na- 
ture of the degrees of freedom in each of these sectors, 
and we start counting them. In terms of spin-1/2 degrees 
of freedom, the original system has 3N spins, the plaque- 
tte operators remove 2F — 2 of them and the nontrivial 
string operators of a given color 2 — x more. These leaves 
N' = 2N — x effective spins in a given sector. Since we 
are not interested in global degrees of freedom, let us 
suppose that the topology is trivial, so that 2 — x = 0. 
Then the resulting value of N' can be easily understood: 
we are only left with hardcore boson degrees of freedom 
(hence the 2N), subject to the constraint of color charge 
preservation (hence the —2). In fact, in each sector we 
can choose a basis with elements labeled by the state of 
the hardcore bosons, so that these become the relevant 
degrees of freedom. As we will see, the dynamics of the 
system transform the hard-core bosons into three families 
of fermions with anyonic statistics between them. More- 
over, we will see that these fermions interact strongly. For 
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arbitrary couplings this only provides us with a picture to 
understand the system, since the emergent quasiparticles 
could be very different. However, in the regime analyzed 
in section fVTl the anyonic fermions come to life explicitly. 

In order to check the statistics of the effective quasi- 
particles emerging from hardcore boson degrees of free- 
dom it suffices [12] to study the hopping terms in ([13]) . 
Strictly speaking, the situation here is not the same as 
in [12], as we are not studying quasiparticles. More close 
is the approach in [l[ , where the emphasis is done on the 
properties of string operators. Now, we first note that 
the hopping terms t c c and t c c that appear in (|13p are only 
enough to hop a c-boson around a c-plaquette. We need 
also composite hoppings: t c c = uiu c ^ c — u^u=' c hops a 
c-boson from a c-plaquette to another. Notice that this 
notation completely agrees with (fT5|) . Consider a state 
with only two c-boson excitations, located at two sites 
separated respectively by a c and a c link from a given 
reference site, as in Fig. [2tc). We may then consider a 
hopping process, indicated with numbers in the figure, in 
which the c-bosons are exchanged in such a way that lo- 
cal contributions cancel [l[ . The net effect of exchanging 
the c-bosons is 



l c b c,c b c l c,c l c b c,c 



, ^.y JU ^zzxx\2 

< r r ,e T T : c T T ,c) 



(27) 



showing that the hardcore-bosons give rise to 
fermions [l^. 

These fermions carry a nontrivial color charge, and 
thus we have three families of them. They are not 
free but interacting, as follows from the existence of the 
u c c terms in the Hamiltonian, which correspond to a 3- 
fcrmion interaction vertex, see Fig. [2fb). In fact, the 
existence of this vertex indicates that we are not dealing 
just with fermions, but rather with anyons: three conven- 
tional fermions cannot fuse into the vacuum. In addition, 
fermions must be coupled to a nontrivial gauge field [13] . 

In order to understand these two issues, let us con- 
sider a process in which a c-fermion is carried around a 
region A, as in Fig[2^d). For clarity, we assume that no 
fermions but that to be transported are present along the 
boundary of A. The hopping process yields a phase 



<t> c x = xc(gx) 



(28) 



where denotes the number of c-fermions inside A and 
the number of 2-colex plaquettes inside A with a num- 
ber of edges that is a multiple of four. We conclude that 
each family of fermions carries a different representation 
of a Z2 x Z2 gauge group with values g\ dictated by the 
vortex states. In addition, fermions of different color have 
semionic mutual statistics: they pick a -1 phase when one 
of them winds around the other. Thus, fermions not only 
interact via virtual fermion exchanges, but also topolog- 
ically. 



VI. EFFECTIVE TOPOLOGICAL COLOR CODE 

We now turn to a perturbative study of the regime 
J z > 0, \J X \, \Jy\ •C J z , for which the previous bosonic 
mapping is specially suited. As in [13j |. we apply the 
PCUTs method [1J| (perturbative continuous unitary 
transformations). This method produces an effective 
Hamiltonian H c g, at a given perturbation order, such 
that [H e g, Q] = 0. We are specially interested in the low- 
energy physics, that is, the Q = sector of the Hamil- 
tonian (notice that the sectors here have nothing to do 
with the ones discussed in the previous section!). For 
each c-plaquette n let us set B% := B°, B% := B%. Then 
since Q — we have 



(29) 



For the particular lattice of Fig. [TJ the 9th order pertur- 
bative calculation yields, up to a constant, the effective 
Hamiltonian 



with 



h — — I J J 1 3 
z q ,J x u y I 



0(J 7 ), 
55489, 



i^x _ l^y _ uwju I j j ,3 

JJ 3 _ 13824 y ■ 



(30) 

(31) 
(32) 



The Hamiltonian (15T))) describes a color code model on 
the effective spins [6|. Its ground state is the vortex free 
sector, where Xc{g-n) = 1 f° r an Y plaquette n. Excita- 
tions are vortices, gapped and localized at plaquettes 
with B% = — 1 and/or J5^ = — 1. At higher orders in 
perturbation theory the terms of the effective Hamilto- 
nian take the form of products of vortex operators, which 
gives rise to vortex interactions 

Color code models have a topological order described 
by a Z2 x Z2 quantum double[15|]. Thus there exist 16 
topological charges in the model, labeled with pairs (<?, x) 
with g G Z 2 x Z 2 and \ an irrep of this group. The 
trivial charge or vacuum sector is (e, Xe)i ^ ne charges 
(c,Xe), (e,Xc), (c,Xc) are bosons and (c, Xe), (c,Xs) are 
fermions. As for the mutual statistics, moving a (g,x) 
charge around a (g', x') charge gives a topological phase 
x(.g')x'(ff)- Regarding fusion rules, a (g,x) charge and a 
(d't xO charge form together a (gg', xx') charge. Notice 
how the three fermions (c, Xe) form a family closed under 
fusion. The same is true for the three (c,xs) fermions. 
Two fermions from different families have no topological 
interaction. This property will turn out to be important 
in the understanding of the model. 

There are several possible ways to label the excitations 
of the gapped phase ([30)1 . Our choice is that an excitation 
at a c-plaquette tt has (c, Xe) charge if — B% = B% = 1, 
(e, Xc) charge if B% = —B% = 1 and (c,Xc) charge if 
Bl=Bt = -\. 
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A. Beyond the low-energy sector 

We now turn our attention to the high-energy quasi- 
particles, that is, those which contribute to the quasi- 
particle counter Q. From the analysis of section [V] we 
now that these are three families of interacting anyonic 
fcrmions. The virtual processes of these fermions create 
the gap for low-energy excitations, as well as interactions 
between them. 

High-energy quasiparticles belong to a particular topo- 
logical charge superselection sector. This correspondence 
between low- and high-energy excitations follows from 
the global constraints (12"2")) . (f2U)) . which have local con- 
sequences. Namely, the local creation of a single o 
fermion must be accompanied by the change of sign of 
several vortex operators, with the net result of the cre- 
ation/annihilation of a low-energy (c, xb) charge. This 
topological charge has exactly the properties needed to 
agree with the fermion fusion rules and the results (|27|) 
and (|28j) derived before. 

Let us now take a closer look to the sector with Q = 1 , 
which involves no interactions. Using the PCUTs method 
one can derive the effective Hamiltonian for the gapped 
phase in the sector with a single c-fermion, for a given 
color c. At first order the c-fermion just hops around 
c-plaquettes, so that for J = J x — J y we get a — 2J 
contribution to the energy gap coming from this orbital 
motion. Going to second order we get a non-flat disper- 
sion relation that corresponds to the triangular lattice 
formed by c-plaquettes. The gap, at this order, is given 
by 1 - 2J - J 2 /2 and thus it closes at J ~ 0.45. This 
is just an approximate estimation, since we are omitting 
all fermion interactions and, perhaps more importantly, 
we are taking J ~ J z . However, it is to be expected 
that as the couplings J x ~ J y grow in magnitude the gap 
for high-energy fermions will reduce, producing a phase 
transition when the gap closes. Such a phase transition 
resembles the anyon condensations discussed in [l6[ , [TtJ • 
Here we are dealing with a condensation of three anyonic 
fermions that does not fit those examples, but part of the 
physical picture could be similar. There exist three non- 
trivial topological charges that do not interact with the 
three condensed anyons. We would thus expect a residual 
topological order in the new phase, with three nontriv- 
ial topological charges with relative semionic statistics. 
These charges would then be responsible for the 4 s ex- 
act topological degeneracy. This conjecture is to some 
extent supported by the validity in all phases of closed 
string integrals of motion which, as we will see next, rep- 
resent processes related to the three topological charges 
invisible to the moving fermions. 

B. Invisible topological charges 

It is natural to expect that open c-string operators 
will create, transport and destroy a particular topologi- 
cal charge among the 16 possible ones. But since string 



operators commute with the hopping terms in (|13p . the 
corresponding three charges must have bosonic relative 
statistic with respect to the moving high-energy fermions. 
That is, they must be 'invisible' for them. This is indeed 
the case, since, as can be easily checked, c-strings move 
(c, Xe) charges. These are three fermionic charges with 
semionic mutual statistics. 

Such an invisibility feature for several charges with 
nontrivial mutual statistics is not present in the Kitaev 
model and has a potential qualitative advantage from an 
experimental perspective. In this regard, in Q it was dis- 
cussed how the high-energy fermions of the honeycomb 
model can spoil attempts to demonstrate the topologi- 
cal order. The fact that in the present model there ex- 
ist topological charges invisible to high-energy fermions 
could simplify the operations needed to show the appear- 
ance of topological phases, because a process that only 
involves the invisible charges is nothing but a product of 
string constants of motion. 



VII. CONCLUSIONS AND PROSPECTS 

We have presented a quantum lattice Hamiltonian in 
2D spatial dimensions that models many relevant and in- 
teresting features regarding systems with quantum topo- 
logical properties. An instance of this is the existence 
of string-nets configurations which are constants of mo- 
tion of the Hamiltonian and are related to fermionic pro- 
cesses. This is a non-perturbative result, thus valid for 
any regime of coupling constants in the model. 

Another prominent result is the emergence of fermionic 
quasiparticles with non-trivial relative statistic, that is, 
anyonic fcrminons. These fcrmions are not free but 
highly interacting, which makes difficult an analytical ap- 
proach. 

The model we have introduced deserves further study 
in many directions that we have motivated throughout 
our work and they are beyond the scope of this paper. 
We can mention hereby several prospects for future work. 

The Hamiltonian Jl} has the symmetries of the lattice 
in Fig. []Ja). We could consider a more general Hamil- 
tonian, in such a way that some of these symmetries are 
lost and we are left only with those of Fig. QJb). This 
amounts to explicitly break the color symmetry of the 
model by substituting the three couplings J w with nine 
couplings J£,. The couplings correspond to the blue 
bonds that connect a c and a c vertex, and the couplings 
J x and Jy correspond respectively to the red and green 
bonds lying on c-plaquettes, see Fig.[IJb). This gives rise 
to a richer phase diagram and physics. For example, in 
the effective color code phase the differences among the 
J° couplings amounts to different gaps for r-, g- and b- 
fcrmions. Then different phase transitions will produce 
as each of the gap closes. 

We have made a clear connection of the model with 
the topological color code, which appears as a gapped 
phase. We have found perturbative indications for the 
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existence of other phases. Thus, it is desirable to con- 
tinue the study of the whole phase diagram of the model 
using other types of techniques, including numerical sim- 
ulations. 

Another aspect that deserves further study is the prop- 
erties of these models for quantum computation in a more 
explicit way, since here we have only focused on their 



properties as far as topological order is concerned and its 
connections with topological color codes. 
Acknowledgements We thank X.-G. Wen and J. Vidal for 
useful discussions and correspondence, respectively. We 
acknowledge financial support from a PFI grant of EJ- 
GV, DGS grants under contract FIS2006-04885 and the 
ESF INSTANS 2005-10. 



[1] A. Yu. Kitaev; Ann. of Phys. 321, 2-111 (2006). 
[2] X.-Y. Feng,G.-M. Zhang, T. Xiang, Phys. Rev. Lett. 98, 
087204 (2007). 

[3] H. Yao, S.A. Kivelson; Phys. Rev. Lett. 99, 247203 (2007). 
[4] S. Yang, D.L. Zhou, CP. Sun; Phys. Rev. B76, 180404 
(R) (2007). 

[5] K.P. Schmidt, S. Dusuel, J. Vidal; Phys. Rev. Lett. 100, 
177204 (2008). 

[6] H. Bombin, M.A. Martin-Delgado; Phys. Rev. Lett. 97, 
180501 (2006). 

[7] M. Kargarian, H. Bombin, M.A. Martin-Delgado; in 
preparation. 

[8] A. Micheli, G.K. Brennen, P. Zoller; Nat. Phys., 2, 341 
(2006). 

[9] V. Lahtinen, Graham Kells, Angelo Carollo, Tim Stitt, 



Jiri Vala, Jiannis K. Pachos Ann. of Phys. 323, 2286 
(2008). 

[10] H. Bombin, M.A. Martin-Delgado; Phys. Rev. B 75, 
075103 (2007). 

[11] M. Oshikawa, T. Senthil; Phys. Rev. Lett. 96, 060601 (2006). 
[12] M. Levin, X.-G. Wen; Phys. Rev. B 67, 245316 (2003). 
[13] J. Vidal, K.P. Schmidt, S. Dusuel: larXiv:0809.1553'y i. 
[14] C. Knetter, G. Uhrig; Eur. Phys. J. B13, 209 (2000) 
[15] F.A. Bais, Peter van Driel, Mark de Wild Propitius 

Phys.Lett. B280, 63 (1992). 
[16] H. Bombin, M.A. Martin-Delgado; Phys. Rev. B 78, 

115421 (2008). 

[17] F.A. Bais, J.K. Slingerland; larXiv :0808.0627 



